Signal-to-noise ratio of Gaussian-state ghost imaging 
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The signal-to-noise ratios (SNRs) of three Gaussian-state ghost imaging configurations — 
distinguished by the nature of their light sources — are derived. Two use classical-state light, 
specifically a joint signal-reference field state that has either the maximum phase-insensitive or 
the maximum phase-sensitive cross correlation consistent with having a proper P representation. 
The third uses nonclassical light, in particular an entangled signal-reference field state with the 
maximum phase-sensitive cross correlation permitted by quantum mechanics. Analytic SNR ex- 
pressions are developed for the near-field and far-field regimes, within which simple asymptotic 
approximations are presented for low-brightness and high-brightness sources. A high-brightness 
thermal-state (classical phase-insensitive state) source will typically achieve a higher SNR than a 
biphoton-state (low-brightness, low-fiux limit of the entangled-state) source, when all other system 
parameters are equal for the two systems. With high efficiency photon-number resolving detectors, 
a low-brightness, high-flux entangled-state source may achieve a higher SNR than that obtained 
with a high-brightness thermal-state source. 

PACS numbers; 42.30.Va, 42.50.Ar, 42.50.Dv 



I. INTRODUCTION 

Ghost imaging is a transverse imaging modality that 
exploits the cross correlation between two photocurrents, 
arising from the detection of two distinct but highly- 
correlated optical beams, to image an object [ll,[3,ii,[|. 
One beam illuminates the object prior to detection by 
a single-pixel (bucket) detector, while the other under- 
goes only free-space diffraction before being detected 
by a high spatial resolution (scanning pinhole or CCD 
camera) detector. Ghost imaging was initially demon- 
strated with biphoton-state light obtained from spon- 
taneous parametric downconversion which requires 
a quantum-mechanical description for its photodetec- 
tion statistics. Subsequent demonstrations with thermal- 
state light [3, [1] , which admits to a semiclassical inter- 
pretation of its photodetection statistics, has generated 
interest in applying ghost imaging to remote-sensing ap- 
plications 0. 

In [lo| we developed a Gaussian-state framework for 
the analysis of ghost imaging that provides a unified 
treatment of biphoton-state and thermal-state illumi- 
nation There we also introduced a classical- 
state source of maximum phase-sensitive cross correlation 
whose ghost-imaging characteristics are most similar to 
those obtained with biphoton-state illumination. For all 
of these sources we determined the near-field and far- 
field image resolution they afford in lensless ghost imag- 
ing, and we quantified the low cross-correlation contrast 
seen with classical Gaussian-state sources and the signif- 
icant advantage, in this regard, that accrues from use of 
the biphoton state. We did not, however, address the 



*Prcsent address: Jet Propulsion Laboratory, Pasadena, California 
91109, USA 

t Electronic address: [baris. i.erkmen@jpl.Dasa.gov| 



signal-to-noise (SNR) behavior of these ghost imagers, 
although we noted the relevance of having high cross- 
correlation contrast to achieving high SNR. This relation 
between SNR and contrast (visibility) has been noted in 
other earlier treatments of ghost imaging as well, e.g., 

a. 

The low cross-correlation contrast of classical-state 
ghost images — which originates from the appreciable 
featureless background in which the desired image 
is embedded — is easily remedied by forming cross- 
covariance images, rather than cross-correlation images 
d, d, This can be accomplished by AC-coupling 

the photocurrents into a correlator, as was done in [isj . 
or by background subtraction. Nevertheless, these tech- 
niques do not eliminate the noise (shot noise and excess 
noise) associated with the featureless background, which 
affects the integration time needed to obtain an accu- 
rate cross-covariance estimate. Therefore it is important 
to quantify the performance of classical and quantum 
ghost imagers via their signal-to-noise ratios. Further- 
more, pursuing closed-form analytic expressions for their 
SNRs is beneficial in identifying the most critical source 
and detector parameters that impact image quality. Sev- 
eral valuable contributions have been made toward this 
end [U, [IBl , but the complexity of the variance expres- 
sion for the image estimate has thus far prevented a rig- 
orous treatment of ghost-image SNR behavior. 

In this paper, we shall utilize our previously-developed 
Gaussian-state framework to derive tractable analytical 
expressions for the SNRs of three lensless ghost imagers, 
whose configurations are distinguished by the nature of 
their light sources. Two use classical-state light, specif- 
ically a joint signal-reference field state that has either 
the maximum phase-insensitive or the maximum phase- 
sensitive cross correlation consistent with having a proper 
P representation. The third uses nonclassical light, in 
particular an entangled signal-reference field state with 
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FIG. 1; (Color online) Lensless ghost imaging configuration. 



the maximum phase-sensitive cross correlation permit- 
ted by quantum mechanics. Because the low-flux, low- 
brightness limit of the last state reduces to vacuum plus 
a weak biphoton component, our analysis encompasses 
biphoton-state ghost imagers. The rest of the paper is 
organized as follows. In Sec. [Til we establish our notation 
and list the general assumptions used in our analysis. 
Then, for each source, we develop its ghost-image SNR 
expression and its low-brightness and high-brightness 
asymptotic behavior in both the near-field and far-field 
regimes. In Sec. IIIII we compare the required averaging 
times for each source to achieve a desired SNR value. We 
conclude, in Sec. IIV| with a summary and discussion of 
our results. 



II. ANALYSIS 



for X = S, R, but vanishing phase-sensitive autocorrela- 
tions. Each source will therefore be distinguished by its 
(phase-insensitive and phase-sensitive) cross-correlation 
functions, which will be specified in the subsections that 
follow. In ([3]): P denotes the photon flux of the sig- 
nal and reference; aq is their beam radius; po is their 
coherence radius, which is assumed to satisfy the low- 
coherence condition po <^ oq; and Tq is their coherence 
time. 

The source-plane signal and reference field operators 
both undergo quasimonochromatic, paraxial diffraction 
over L-m-long free-space paths, yielding detection-plane 
field operators Ei{p,t) and E2{p,t) respectively. These 
detection-plane field operators are also in a zero-mean 
jointly-Gaussian state, whose second-order correlation 
functions can be derived from their source-plane counter- 
parts and the free-space Green's function [l^, [l3| . The 
first field, Ei{p,t), illuminates a quantum-limited pin- 
hole photodetector whose photosensitive region is cen- 
tered at the transverse coordinate Pi , and whose photo- 
sensitive area Ai is smaller than the coherence area of 
the impinging field state. The second field, E2{p,t), il- 
luminates a field-transmission mask T{p) that is located 
immediately in front of a quantum-limited bucket pho- 
todetector which collects all light transmitted through 
the mask fisj . 

The photodetectors are assumed to have identical sub- 
unity quantum efficiency rj and no dark current or elec- 
tronic thermal noise. Furthermore, finite-bandwidth 
postdetection filters AC-couple the photocurrents into the 
correlator block from Fig. [U so that a background-free 
ghost image is obtained. For analytic simplicity, we as- 
sume that the composite baseband frequency response of 
the photodetectors and their AC-coupling is given by the 
difference of two Gaussian functions [191 , 



The lensless ghost imaging configuration that we will 
consider is shown in Fig. [T] [l^. Here, Es{p,t)e~^'^°^ 
and ER{p,t)e~^'^°* are scalar, positive-frequency, parax- 
ial, signal (S) and reference (R) source-field operators, 
each with center frequency ujq and ^/photons/m^s units. 
These operators obey the canonical commutation rela- 
tions [17| 



[Em{Pi,ti),EeiP2,t2) 
[Em{Pl,h),El{p2,t2) 



(1) 

SnijS{Pi- P2)S{ti-t2) , (2) 



for m,£ = S, R. The joint signal-reference source states 
we shall employ will all be zero-mean jointly-Gaussian 
states in which the signal and reference beams have iden- 
tical coherence-separable Gaussian-Schell model phase- 
insensitive autocorrelation functions given by 



{Elip„h)E,ip^,t2)) 
2P 



(3) 



(4) 



where VIb is the baseband bandwidth of the detector 
(taken at the attenuation level of the frequency re- 
sponse), r^AT ^ VIb is the stopband bandwidth of the 
AC-coupfing notch around SI = 0, and T[hB{t)] denotes 
the Fourier transform of the composite filter's impulse 
response, hB(t). In order to minimize suppression of the 
baseband photocurrent fluctuations — whose cross corre- 
lation yields the ghost image — the notch bandwidth will 
be taken to be much smaller than the bandwidth of the 
impinging fields, i.e., ^nTq <C 1 will be assumed in all 
that follows. 

The ghost image at the transverse coordinate Pi is 
formed by time-averaging the product of the detector 
photocurrents, which is equivalent to a measurement of 
the quantum operator 



CiPi) = ^ [ ^ dth{t)h{t) 



(5) 
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where [2C 



?mW=9y"dT^ dpEi{p,T)E'^{p,T)hB{t-T), (6) 

for m = 1,2, with q being the electron charge. Am the 
photosensitive region of detector m, and T/ the duration 
of the averaging interval. The field operators appearing 
in these photocurrent operators are 



(7) 



[ ^T{p)E2{p,t) + VI -?7|T(p)pi?vac.(p,t), 

for m = 1, 2, where the {i?vac„ (p, i)} — which are needed 
to ensure commutator preservation — are in their vacuum 
states. 

The C{pi) measurement yields an unbiased estimate 
of the ensemble-average equal-time photocurrent cross- 
correlation function 

{C{p,))^{h{t)h{t))^q'rj'A^ 

dp Jdui j du2 hsit - ui)hB{t - U2)\T{p)\^ 

X {E\{p^,ui)El{p,U2)Ei{pi,ui)E2{p,U2)) , (8) 

where we have approximated the integral over the pin- 
hole detector's photosensitive region as the value of the 
integrand at pi times the photosensitive area Ai . Evalu- 
ation of {C{pi)), for the Gaussian-state sources we shall 
consider, can be accomplished along the lines established 
in To find the ghost image signal-to- noise ratio at 
the point Pi, 



SNR 



i^cHp,))' 



(9) 



where AC{p) = C{p) — {C{p)), all that remains is to 
evaluate the variance term appearing in the denominator. 
We have that the variance term obeys 

(AC2(Pi)> = 

— 2 / dt j du {ii{t)i2(t)ii{u)i2iu)) 

J-Ti/2 J-Ti/2 



(10) 



This expression reveals the primary challenge in evaluat- 
ing the measurement variance: the fourth-moment of the 
photocurrents in the integrand is an eighth-order moment 
of the field operators. Fortunately, the moment-factoring 
theorem for Gaussian-state optical fields [2l|, [l^l — which 
we used in [IQ] to find {C{pi)) — allows all field moments 



to be expressed in terms of second-order moments. Be- 
cause this procedure is straightforward but tedious, we 
shall confine our discussion here to a detailed descrip- 
tion of the simplification procedure, rather than a lengthy 
derivation. 

First, we express the integrand on the right-hand side 
of (jlOp in terms of the field-operator moments, as we 
have done in ^ for the mean. We then use the com- 
mutator relations H]), ([2]) to put the integrand into 
normal order. This procedure yields the sum of four 
normally-ordered moments: one eighth-order moment, 
two sixth-order moments, and one fourth-order moment. 
Next, the Gaussian-state moment- factoring theorem is 
applied to each term, replacing higher-order moments 
with expressions that depend only on the second-order 
moments of the fields. Note that the nonzero terms in 
the moment-factored expression depend on whether the 
source of interest has nonzero phase-sensitive or phase- 
insensitive cross-correlation functions. Finally, employ- 
ing the coherence-separability of the correlation func- 
tions, the spatial and temporal integrals in each term are 
evaluated separately. It is relevant to note that many 
temporal integrals vanish due to our use of AC-coupling, 
i.e., because Hb{0) = 0. Moreover, the symmetry prop- 
erties of correlation functions can be used to group some 
nonzero terms so that the final variance expression is a 
sum of only eight terms. We now proceed with the details 
for each of the three sources under consideration. 



A. Thermal-State Light 

Lensless ghost imaging with thermal-state light usu- 
ally derives its signal and reference sources from 50-50 
beam splitting of a single zero-mean Gaussian-state beam 
possessing a phase-insensitive autocorrelation function 
but no phase-sensitive autocorrelation function, see, e.g., 
[l^ . Taking the post-splitter signal and reference fields 
to have the Gaussian-Schell model autocorrelations from 
([3]) , it follows that these fields have the maximum phase- 
insensitive cross correlation, given by 



{El{p„ti)ER{p^,t2)) = 



2P 



e-(IPil'+IP2l')/a^-|P2-PilV2p?e~(*^"*i)'/2T,?, (11) 



and a vanishing phase-sensitive cross correlation, viz., 
{Es{p„h)ER{p^,t2)) ^0. 

Let us begin our thermal-state SNR analysis with near- 
field operation, wherein kopoao/2L 3> 1 prevails, with 
fco = loq/c being the wave number associated with the 
center frequency loo. In this regime, the detection-plane 
correlation functions are approximately equal to those of 
the source. With all auto- and cross-correlation functions 
specified, evaluating the spatial and temporal integrals 
in the moment-factored variance expression is a straight- 
forward exercise. For the spatial integrals, we assume 
that flo exceeds the transverse extent of the transmission 
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mask by an amount sufficient to permit the approxima- 
tion e~l''l ^°'o\T{p)\ « |r(p)|. For convenience, we define 



dp\T{p) 



{12) 



which we will regard as the effective area of the trans- 
mission mask. Our A'j. interpretation follows by anal- 
ogy with the case of a binary (|r(p)| G {0,1}) mask, 
for which is the area over which — 1. With 

this interpretation we have that Nrpj p\ is the number of 
spatial resolution cells in the ghost image We also 
note that the small-pinhole approximation introduced in 
the previous section requires p\l A\ ^ 1 for its validity in 
near-field operation. Finally, we identify the two assump- 
tions employed in evaluating the variance expression's 
temporal integrals: T/ ^ Tq and fisT/ ^ 1. Neither of 
these averaging-time conditions is at all surprising. The 
former states that we must average over many source co- 
herence times to form a high-quality ghost image. The 
latter states that we must average over many photode- 
tector response times to achieve this same purpose. 

Within the near-field regime — and subject to the as- 
sumptions given in the preceding paragraph — we will 
evaluate the ghost-imaging SNR behavior that prevails 
under narrowband and broadband illumination condi- 
tions. A source state is said to be narrowband if Sis To ^ 
1, so that the coherence time of the impinging field state 
To greatly exceeds the ^l/fis integration time of the 
photodetectors. Conversely, a broadband source state is 
one that satisfies fisTo ^ 1, so that the source's co- 
herence time is much shorter than the photodetector's 
integration time. 

For a narrowband source and near-field ghost imaging 
with thermal-state light when A'rp/ ^ 30 (e.g., the 2D 
image consists of 10 x 10 or more resolution cells), we 
find that the signal-to-noise ratio is 



SNR = 



\T{p,)\^Ti/To 



inpi)P 



4:7rpf,\Tip,)\^ 



v/¥OBTopg|T(pi)P 
16V2Ai?72l2 



,(13) 



where X = PToPq/uq is the source brightness, i.e., the 
source's average number of photons per spatiotemporal 
mode. As expected, this SNR expression grows linearly 
with increasing averaging time Tj, behavior that will be 
seen in all the cases we will consider in this paper. More 
importantly, we can give physical interpretations to the 
terms in its noise denominator that dominate in low- 
brightness and high-brightness operation. 

We have chosen to use quantum photodetection theory 
to derive all the SNR expressions in this paper. How- 
ever, as shown in [10], quantitatively identical formu- 
las follow from semiclassical photodetection theory when 
the signal- reference state is classical, i.e., when it has 
a proper P representation. The thermal state is classi- 
cal. It is thus appropriate to replace the photocurrent 



operators {im{t)} with classical photocurrents {im{t)} 
that, owing to the assumed AC-coupling, are zero-mean 
random processes comprised of a shot-noise component, 
arising from the discreteness of the electron charge, plus 
an excess noise component, which is proportional to the 
fiuctuations in the photon fiux illuminating the detec- 
tor. The variance contributions generated by these pho- 
tocurrent components scale differently with increasing 
source brightness. As a result, we can identify the left- 
most and right-most terms in the noise denominator of 
(fT5)l — which are the noise terms that dominate at high 
and low source brightness, respectively — as being nor- 
malized variance contributions coming from excess noise 
alone and from shot noise alone, while the middle terms 
arise from beating between excess noise and shot noise. 
Thus, as the source brightness grows without bound, the 
SNR from ([13]) increases until it saturates at its maxi- 
mum value, 



SNR„,ax = 



(14) 



which is limited by the excess-noise term [23|. Roughly 
speaking, SNR^ax equals the number of source coherence 
times in the averaging interval divided by the number of 
spatial resolution cells in the image and multiplied by the 
square of the object's intensity transmission. Note that 
Pq/A'j, is the image contrast for DC-coupled ghost-image 
formation in the near field with narrowband thermal- 
state light [lOj. Hence, the SNR of AC-coupled, high- 
brightness, thermal-state ghost imaging is proportional 
to the image contrast realized using the same setup with 
DC-coupling. 

At very low source brightness the ghost-image SNR 
obtained with a narrowband thermal-state source will be 
controlled by the shot-noise contribution to its noise de- 
nominator. In this case (fT^ reduces to 



16V2 TiT^PA, 2 
SNR = — - — jT]I\T{p^)\ . 

y/n Jo "sao 



(15) 



In Fig. 2(a) we have plotted the narrowband thermal- 
state ghost imaging SNR from — along with its high- 
brightness and low-brightness asymptotes from (|14p and 
psp — for several narrowband ghost imaging scenarios. 

We now turn our attention to broadband sources, 
which satisfy fisTo < 1. When A'rp/pl > 12 f^sTo holds, 
we find that 



SNR = 



\Tip,)\^Tj/To 



2|T(pi)P 



87rpg|T(pi)|4 



,(16) 



where, once again, the left-most and right-most terms in 
the noise denominator are due to excess noise alone and 
shot noise alone. Here too SNR increases with increasing 
source brightness until it reaches its maximum value. 
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nBTi^\T{p,)\\ 



(17) 
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FIG. 2: (Color online) Thermal-state ghost imaging SNR, 
normalized by Ti /To, plotted as a function of source bright- 
ness X = PTopl/al, for \T{p^)\ = 1, A'^/pl = 10^ pi/ A, = 
10 and ?7 = 0.9. Various fisTo values are shown in the (a) nar- 
rowband and (b) broadband limits. Dash-dotted lines repre- 
sent low-brightness asymptotes, and dashed lines correspond 
to high-brightness asymptotes. 



where it is limited by excess noise alone. Aside from 
insignificant numerical factors, the broadband SNRmax 
differs from the narrowband SNRmax only through re- 
placement of 1 /To from the narrowband expression with 
fis in the broadband expression. This replacement is to 
be expected. In the narrowband case T/ /Tq is the number 
of photocurrent coherence times that are being averaged 
by the correlator. This is because the narrowband condi- 
tion r^sTo 3> 1 ensures that the photon- flux fluctuations 
are not affected by the photodetector's baseband band- 
width limit. However, under the broadband condition, 



^bTq <C 1, the photocurrent fluctuations have a much 
longer (^l/rig) coherence time than that of the photon 
flux illuminating the detectors, so it is f^sT/ that appears 
in the broadband SNRmax formula. 

At very low source brightness, the SNR of the broad- 
band thermal-state ghost imager becomes limited by shot 
noise alone and is given by 



SNR = 



(18) 



Figure 2(b) shows several plots of broadband thermal- 
state ghost imaging SNR, together with its high- 
brightness and low-brightness asymptotes. 

Thus far we have concentrated on the near-field SNR 
behavior with a thermal-state source. Our results, how- 
ever, are easily converted to the far-field regime, in which 
koPoiof^L <^ 1. In order to obtain the far-field SNR we 
must first propagate the second-order correlation func- 
tions from the source plane to the detection planes. For 
Gaussian-Schell model correlation functions this trans- 
formation is a simple replacement of ao by = 2L/kopo, 
and Pa hy PL = '^L/koao, when quadratic phase fac- 
tors that do not affect ghost-image formation are omitted 
[2^ . Therefore, all of our near-field thermal-state SNR 
results can be converted to corresponding far-field results 
by making these parameter value changes. 



B. Classically-correlated phase-sensitive light 

Let us now consider a source state that has the maxi- 
mum phase-sensitive cross correlation permitted by clas- 
sical physics, given the auto-correlation functions in ([3]), 
but has no phase- insensitive cross-correlation, i.e.. 



2P 



-(|Pil' + |P2l')/io-|P2-PilV2Pop-(t2-ti)V2T= 



(19) 



and {Eg{pi,ti)Eji{p2,t2)) = 0, where we have assumed 
the phase-sensitive cross correlation function is real val- 
ued i2|. 

In the near-field regime, which is now given by 
kopQ/2L ^ 1, we can follow the same assumptions 
stated in the previous section for thermal states and ar- 
rive at the same near-field SNR expressions, i.e., (jl3[) 
and (|16p apply to narrowband and broadband near-field 
operation with this classically-correlated phase-sensitive 
source state. Furthermore, the high-brightness and low- 
brightness asymptotes are as given by (fT4|) and ((T5)) . re- 
spectively, for the narrowband case, and by (|17p and 
iHH]), respectively, for the broadband case. Moreover, 
the physical interpretations we provided for the thermal- 
state results continue to apply, without modification, for 
the classically-correlated phase-sensitive source. 

For phase-sensitive coherence, the far-field regime cor- 
responds to ki^a^/2L ^ 1. As stated in the previ- 
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ous section, the detection-plane (phase-insensitive) auto- 
correlation functions are found by substituting for 
ao and pL for po- Propagating the source-plane phase- 
sensitive cross correlation given in (jl9p involves the same 
substitutions, but in addition requires replacing |/32~PiP 



by |P2 



because the far-field ghost image formed 



with phase-sensitive light is inverted. It follows that the 
the far-field SNR expressions are derived from the near- 
field SNR expressions by replacing the source-plane co- 
herence radius po and beam radius Oq by their detection- 
plane counterparts, p^ and respectively, and using 
\T{~p,)\ in lieu of |T(pi)|. 



C. Maximally-entangled phase-sensitive light 

We continue to consider signal and reference beams 
in a zero-mean jointly-Gaussian state with no phase- 
insensitive cross correlation, but now we take the phase- 
sensitive cross correlation to be the maximum permit- 
ted by quantum physics. Because quantum ghost imag- 
ing experiments have used the signal and idler outputs 
from spontaneous parametric downconversion (SPDC) as 
the two source fields, we shall focus on this case here. 
The output field operators of SPDC can be expressed 

as (Ulillii 



-B,„(p, t) = A{p)£„-i{p, t) + Cm {p, t) 



(20) 



for m = S,R^ where |^(p)| < 1 is an aperture func- 
tion representing the finite transverse extent of the in- 
teraction medium, and (p, t) are auxiliary vacuum- 
state operators, so that the Em{p,t) satisfy the free- 
space field commutator relations. The operator-valued 
Fourier transforms of {£rn{p,t))'m = S,R}, denoted by 
{Am{k,V,),m = S,R}, are given by a two-field Bogoli- 
ubov transformation of vacuum-state input field opera- 
tors, CLmik, fi), i.e., 

As{k,n) - 

p{k, n)asik, n) + iy{k, 17)a^(-fc, -fl) (21) 

AR{-k,-n) = 

p{k, n)aR{-k, -n) + v{k, n)al{k, fl). (22) 

Here i'{k, ft) E M. and /i(A;, 51) = 1 -f ii^{k, il) are the 
canonical transformation coefficients. In accordance with 



the Gaussian-Schell model treatment introduced earlier, 
we set ii 



i^{k,n) = 2(27r)i/ 



1/4. /^^OPo„-pf,|fc|V4-To^oV4 



and 



^(p) = cxp{-|p|Vao}, 



, (23) 



(24) 



such that the Em{p, t), for m = S, R, are in a zero-mean 
jointly Gaussian state, with phase-insensitive autocorre- 
lation functions given by ([3]), and the maximum permis- 
sible phase-sensitive cross correlation function, 



(^s(pl,^l)^J^(p2,^2)) - — 



-(|Pll' + |P2l')/ao 



X e 



TTOn 



=-|P2-PllV2pL-(t2-tl)V2To' 



.e-|P2-PilVpoe-(*2-*i)Vro' 



.(25) 



All other second-order moments, i.e., the phase-sensitive 
autocorrelation functions and the phase-insensitive cross- 
correlation function, are zero. It is worthwhile to point 
out that when the source brightness X = PTop'^/a^ ^ 1, 
the first term in the square brackets dominates, and ([25]) 
approaches the classical phase-sensitive cross correlation 
given in (jl9p . However, when X <|; 1, the second term is 
much larger than the first, resulting in a much stronger 
phase-sensitive cross correlation than permitted in a clas- 
sical state. If the brightness is lowered to the limit in 
which there is on average much less than one photon 
in the signal and idler beams, the output of the SPDC 
can be approximated as a dominant vacuum component 
plus a we ak p air of entangled photons, viz., the biphoton 
state [13, [Ig. 

To evaluate the ghost-image SNR in the near-field 
regime (fcoPo/2-^ S> 1) we utilize the same approxima- 
tions we have used for classical phase-sensitive light, now 
with the cross-correlation function from (j25p employed 
in lieu of (|19p when integral expressions are explicitly 
evaluated. In the narrowband limit, ^bTq 3> 1, we find 
that 



SNR: 



1 + 1/V2^X Ti/To 



/2^|T(Pi)|Vg 



i(r, 



47rpg 
3Ai7j 



(26) 



16^2^1 |T(pi)|2,,2I2 



/27rl 



This SNR expression captures the full quantum-to- classical transition seen in ghost imaging with maximally- 
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entangled phase-sensitive light (the output fields from 
SPDC), as the source brightness X is increased. When 
the source is bright, i.e., X ^ 1, the first terms in the 
numerator and the denominator of (|26p are dominant, 
yielding the same SNR as that obtained with narrow- 
band, bright classical maximally-correlated (thermal or 
phase-sensitive) light, i.e., 



SNR: 



(27) 



For dim-source (X ^ 1) ghost imaging the second term 
in the numerator of (|26p dominates the first, and when 
l^(Pi)l ^ 1 the last term in the noise denominator is the 
most significant, yielding an SNR that is linear in photon 
flux, 



SNR: 



16r,r;2pAi|T(pi)|^ 



TT Tn 



(28) 



For example, even with the generous parameters values 
PT^pl/al = 0.01, pI/Ai = 10, and Q^bTq = 10, it is nec- 
essary to have 10^ resolution cells before p6)) deviates ap- 
preciably from the linear dependence on the photon flux 
P for |T(p]^)| ~ 1. In this regime the SNR is limited by 
the very low number of photon pairs detected over a de- 
tector integration time. The SNR achieved with narrow- 
band maximally- entangled phase-sensitive light is plot- 
ted in Fig. 3(a) for several flsTo values. The plots verify 



the linear low-brightness regime and the high-brightness 
saturation towards the classical asymptote. However, as 



shown in the plotted curves, the SNR can exceed the 
bright-source asymptote. When this occurs, there is a 
flnite source brightness that yields the maximum SNR, 
and increasing I beyond this threshold will decrease the 
SNR with increasing photon fiux. For a given set of pa- 
rameters this source-brightness threshold can be found 
easily by solving for the roots of a third-order polynomial, 
which yields the critical points of (pS)) . Although closed- 
form solutions for these roots exist, the expressions are 
too cumbersome to pursue further in this paper. 

If the low- brightness condition of the source (I <C 1) 
is augmented with the low-flux condition 



«1, 



(29) 



then the average number of photons per integration time 
impinging on either detector becomes much less than 
unity. It follows that the photodetectors can be replaced 
with non-photon-resolving photodetectors without ap- 
preciable loss in imaging functionality, thereby render- 
ing the Fig. [1] ghost imaging configuration equivalent to 
biphoton-state ghost imaging with coincidence-counting 
circuitry (instead of photocurrent correlation). Thus, 
narrowband biphoton-state ghost imaging is also gov- 
erned by the linear photon-fiux SNR formula from p8|) 
for \T{p,)\ ^ 1. 

Shifting our attention to the broadband (r^sTo ^ 1) 
limit, we arrive at 



SNR: 



1 + 1/20FX Ti/To 



^nBTop^\T{p^)\^ 



Di. 

I 



( 72 I z^pi/i 



(30) 



2'7|T(Pi)|^ 



87,^Ai|T(pi)|^I== 



where 



V8 



When X 3> 1, the first terms in the numerator and de- 
nominator are dominant. Consequently, the SNR ap- 
proaches the SNR of classical ghost imaging with a 
bright (phase-insensitive or phase-sensitive) maximally- 
correlated broadband source, which is given in (fTT)) . 

On the other hand, if the low-brightness condition, 
X <C 1, and the low-flux condition, as given in (|29p. are 
both satisfied, the right-most term in the noise denomi- 
nator becomes dominant, for |r(p]^)| ^ 1, yielding 



IT "BOq 



for the broadband biphoton-state SNR expression. As in 
the narrowband case, the SNR in this regime suffers from 
^32) the paucity of photon pairs detected within the photode- 
tector's ~l/riB integration time, which is compensated 
by averaging the photocurrent product over many tem- 
poral coherence bins, i.e., employing Tj^Ib ^ 1. Fig- 
ure [3(by| shows plots of ([50]) for several values of fisTo. 
It is seen that the linear photon-flux dependence of the 
SNR extends well beyond the low-flux regime. Further- 
more, the SNR achieved in the low-brightness regime 
(i.e., (1221)) can exceed the bright classical source asymp- 
tote given in (|17p . Similar to what we found for narrow- 
band operation, the broadband SNR has a well-defined 
maximum achieved at finite source brightness. Increas- 
ing the brightness beyond this optimal value reduces the 
(32) SNR, which ultimately converges to the SNR attained 
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with classical sources. All three of these SNR regimes 
are clearly identifiable in the Fig. 3(b) SNR plots. For 
example, for QbTq = 10^^, the SNR is linear in photon 
flux imtil X K, 10"'^, its maximum occurs at X s:! 10~^, 
then the SNR decreases with increasing X until at X « 10 
it converges to the bright-source asymptote. 



(a)Narrowband. 




10' 

(b) Broadband. 



10° 10' 



FIG. 3: (Color online) Nonclassical phase-sensitive Gaussian- 
state ghost-imaging SNR, normalized by Ti/To, plotted ver- 
sus source brightness T = PToPo/oq for a near-field config- 
uration {kopo/2L > 1) with |r(pi)| = 1, A'^/p^ = 10'', 
Po/Ai — 10, and rj = 0.9. Various Q,bTo values are shown in 
the (a) narrowband and (b) broadband hmits. Dash-dotted 
lines represent low-brightness asymptotes and dashed lines 
correspond to high-brightness asymptotes. 



In the far field {koaQ/2L ^ 1) the source-plane phase- 
sensitive cross correlation in (j25p must be propagated to 
the detection planes [l^l before the SNR can be evalu- 
ated. For narrowband sources {^IbTq ^ 1) we get 



SNR = 





(l + 1/V8^X^ 








1 { 1 1 47rp2 \ 






\/2¥|T(-pi)|V| ^ 


I\V\T{-Pi)\^' ' 3Air,J 


^ 16y2Ai|T(-pi)|2^2x2 





(33) 



This result simplifies to in the low-flux (biphoton-state) limit, defined by (|29p 

SNR.i^^^^^M ,34, 
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with |T(— ~ 1, showing that narrowband biphoton- 
state ghost imaging SNR is Hmited by the number of 
photon pairs detected within the detectors' integration 
time. On the other hand, when X ^ 1, the SNR con- 
verges to 



Jo 



SNR = V27r;^^|r(-Pi)r, (35) 
which, as explained in the previous subsection, is equal to 



the bright-source SNR asymptote for far-field ghost imag- 
ing with narrowband classical (maximally-correlated) 
phase-sensitive light. Similar to the near-field scenarios, 
the SNR can have a maximum at finite source bright- 
ness. The full behavior of the far-field biphoton SNR in 



the narrowband regime is shown in Fig. 4(a) 



The SNR in the broadband case {QbTo < 1), on the 
other hand, is given by 



SNR: 



(l-|-l/4V^j)V//ro 



(36) 



^nBTopl\Ti~p^)\-^ 



V2 



D2 



+ 16v'Ai\T{~p,)\-^r^ 



where 



V2 Sirpj 



1 



nsTo 37]AA 477|T(-pi)|2y' 
which simplifies in the biphoton limit to 

SNR=W,!^!^MVPr)P^ 



(37) 



(38) 



for |T(— Pi)| ~ 1. Hence, broadband biphoton-state 
ghost-imaging SNR is limited by the number of photon 
pairs detected within the detectors' ~ integration 
time. The X > 1 SNR simplifies to 



SNR=JlnBTj-^\Ti-p,)\\ 



(39) 



which is identical to the bright-source SNR asymp- 
tote for far-field ghost imaging with broadband classi- 
cal (maximally-correlated) phase-sensitive light. The far- 
field biphoton SNR in the broadband regime is shown in 
Fig. |4(b)] for several flsTo values. 



III. IMAGE ACQUISITION TIMES 

All of our SNR expressions are proportional to the 
cross-correlation averaging time Tj. Consequently it is 
meaningful to compare the averaging times required to 
achieve a desired SNR value with different ghost imagers. 
Let us first assume that all parameters except photon 
flux are equal in the configurations of interest. Because 
the classical sources (the thermal state and the classi- 
cal phase-sensitive state) yield identical SNRs, we shall 
use P^'^-' to denote their photon fluxes, reserving 

pig) for 

the photon flux of the maximally- entangled (nonclassical) 
state. Likewise, the averaging time for the classical-state 



ghost imagers to achieve the desired SNR will be denoted 

(c) 

by Tj ' , while that for the entangled-state source will be 

designated Tj'^\ 

Then, in the near field and with narrowband sources 
we obtain [301 



(9) 



Pi 



8V2 r;2p(9)^^ Ai 



\Tip^)\' 



(40) 



where the classical-state sources are assumed to be bright 
enough to achieve the saturation SNR in but the 

quantum source is limited to low brightness so that the 
nonclassical signature of the source prevails. In the 
biphoton-state (low-flux) limit with \T{pi)\ ^ 1, ([iO]) 



implies t}''^ > tI''\ i.e., the cross-correlation averag- 
ing time required for narrowband biphoton-state ghost 
imagers to achieve a desired SNR value is much longer 
than that for bright classical-state ghost imagers to do 
so, given that all other system parameters are equal. 
In the near field using broadband sources we obtain 



(9) 



Pi 



■.(c) 



4^2 r?2p(<;)^^ Ai 



Tip,)? 



(41) 



Thus, when the quantum ghost imager utilizes a low- 
flux (biphoton-state) source and |P(Pi)| ~ 1, we find 
that rj'-* S> Tj'^-' prevails, i.e., once again the cross- 
correlation averaging time required for ghost imaging 
with broadband bright classical-state light is significantly 
shorter than that for biphoton-state ghost imaging, given 
equal system parameters. Nevertheless, in a very high- 
resolution ghost imaging configuration, high illumina- 
tion flux (p(9U^/17Bag > 1) may be achievable with 
low-brightness (P^'^'ToPo/o-o ^ 1) maximally- entangled 
phase-sensitive Gaussian-state light. In this case (PT|) 
implies that the averaging time for the nonclassical-state 
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(b) Broadband. 



FIG. 4: (Color online) Nonclassical phase-sensitive Gaussian- 
state ghost imaging SNR, normalized by Ti /To, plotted versus 
source brightness I = PToPo/oq for a far- field configuration 
(fcoag/2L < 1) with |T(-pi)| = 1, A't/pI = 10^ pg/Ai = 10 
and 77 — 0.9. Various SlsTo values are shown in the (a) nar- 
rowband and (b) broadband limits. Dash-dotted lines repre- 
sent low-brightness asymptotes and dashed lines correspond 
to high-brightness asymptotes. 



ghost imager can be shorter than that for the classical- 
state ghost imager [3l|. For example, |r(pi)| — 1, 
A't/pI = 104, ^bTo = 10-2 p('?)rop2/a2 ^ 10-3^ and 



pI/Ai = 10 will yield t)"^^ w r^'VlOO?]^. 

Finally, we compare ghost imaging with a broadband 
biphoton-state (low-brightness and low-flux) to that with 
a bright narrowband classical state. Denoting the param- 
eters specific to the classical and quantum sources with 



.(c) 



the superscripts (c) and (g), we obtain 

- ^ n^^^a?, pI |T(pi)|2 



B "0 



pI 



(42) 



Because the last factor on the right is typically less than 
unity, for |T(p]^)| 1, whereas the remaining factors are 
greater than unity, the cross-correlation averaging time 
may be favor either source state. As an example, con- 
sider \T{p,)\ - 1, A'^/pI = 10^, P^'^^tI,''^ p^/al = IQ-', 
pI/Ai = 10. The n a biphoton-state source with 1 THz 
bandwidth and a 1 MHz thermal-state source will re- 
sult in rj''^ w 4 X IQ-'^rf^ /jf, which shows that the 
biphoton-state imager enjoys an enormous advantage in 
averaging time as compared to the classical-state imager 
when the quantum efficiency is not unduly low. 

So, ghost imaging with bright classical thermal or 
phase-sensitive states affords a shorter averaging time 
to reach a desired SNR value than does a biphoton- 
state ghost imager, given all of the remaining parame- 
ters governing the two ghost imaging systems are equal 
and |r(pj^)| ~ 1. However, ghost imaging with low- 
brightness, but high-fiux quantum sources may achieve 
shorter averaging times than bright classical sources, if 
high quantum efficiency photon-number resolving pho- 
todetectors are employed. Finally, if ghost imaging with 
a broadband biphoton state is compared to that with a 
narrowband bright classical state, the integration time 
may favor either source, depending on the ratio of the 
achievable source bandwidths, number of resolution cells 
in the image, and the biphoton-state source brightness. 



IV. DISCUSSION 

We have presented a detailed SNR analysis for three 
Gaussian-state ghost imaging configurations. Two used 
classical-state light, specifically a joint signal-reference 
field state that has either the maximum phase- insensitive 
or the maximum phase-sensitive cross correlation con- 
sistent with having a proper P representation. The 
third used nonclassical light, in particular an entangled 
signal-reference field state obtained from SPDC, with the 
maximum phase-sensitive cross correlation permitted by 
quantum mechanics. Our analysis concentrated on the 
narrowband and broadband limits in both the near field 
and the far field. Because the conclusions from our anal- 
ysis apply, in identical ways, to both the near-field and 
far-field regimes, we shall omit references to these regimes 
in what follows. 

We found that classical-state ghost imager SNRs 
saturate — with increasing source brightness — to maxi- 
mum values that are inversely proportional to the num- 
ber of resolution cells on the imaged object. In this 
high-brightness limit the SNR is thus proportional to the 
contrast achieved in DC-coupled operation. The contrast 
can be improved by removing the featureless background 
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via AC-coupled photodetectors or background subtrac- 
tion, but SNR improvements, at high source brightness, 
require increasing the cross-correlation averaging time. 

Biphoton-state ghost imagers were shown to have 
SNRs that are typicaUy proportional to their low pho- 
ton flux. Hence, for such imagers, increasing the pho- 
ton flux of the source is of prime importance. How- 
ever, the SNR gain derived from such increases is not 
unbounded. Our analysis revealed that the SNR re- 
alized with a nonclassical-state source with low bright- 
ness but high flux typically has a well-defined maximum, 
after which increasing flux reduces the SNR. The in- 
verse relation between SNR and photon flux may seem 
counterintuitive, but it is consistent with the fact that 
the SNR must approach the bright classical-state SNR 
as the source brightness increases beyond unity. This 
classical-state limit, however, is in general lower than 
the maximum SNR achieved by the nonclassical-state 
source in the low-brightness regime. Hence, in these 
cases, the SNR achieved with nonclassical phase-sensitive 
light must have a decreasing trend as source brightness 
increases without bound. 

To assess the performance achieved by different sources 
we compared their image acquisition times, i.e., the 
cross-correlation averaging times needed to achieve a 
predefined target value for SNR. We showed that with 
equal bandwidth sources, and all system parameters be- 
ing equal unless otherwise noted, bright classical-state 
ghost imagers typically reach the desired SNR value with 
a much shorter averaging time than that needed by a 
biphoton-state ghost imager. Therefore, although the 
biphoton-state yields images with high contrast even in 



DC-coupled operation, because the biphoton imager is 
photon starved the total time duration it requires to 
accumulate the ghost image far exceeds that necessary 
with a bright classical-state source. Nevertheless, we 
saw that there is a broadband, low-brightness, high-flux 
regime of nonclassical phase-sensitive light that may get 
by with much shorter cross-correlation averaging times 
than those needed by bright classical-state imagers. The 
notable drawback to reaping this quantum advantage, 
however, is the necessity for high quantum efficiency 
photon-number resolving detectors. 

In conclusion, Gaussian-state analysis provides a ro- 
bust and versatile framework for answering some of the 
most fundamental questions associated with developing 
practical ghost imagers for remote sensing applications. 
In this paper, we have used this framework to study 
the SNR behavior of ghost imagers with source states 
that encompass those that have been used in proof-of- 
principle ghost imaging experiments. Our analysis un- 
ambiguously identifies the key parameters that limit SNR 
behavior. For high-brightness classical-state ghost imag- 
ing it is the number of resolution cells in the image, 
whereas for the biphoton-state case it is the low photon 
flux of the source. 
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